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Abstract. We discuss the distribution of radial velocities of galaxies belonging to the 
Local Group. Two independent samples of galaxies as well as several methods of reduction 
from the heliocentric to the galactocentric radial velocities are explored. We applied the 
power spectrum analysis using the Hann function as a weighting method, together with 
the jackknife error estimation. We performed a detailed analysis of this approach. The 
distribution of galaxy redshifts seems to be non-random. An excess of galaxies with radial 
velocities of ~ 24km • s~^ and ^ 36km • s~^ is detected, but the effect is statistically 
weak. Only one peak for radial velocities of ~ 24 km ■ s^^ seems to be confirmed at the 
confidence level of 95%. 
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1 INTRODUCTION 



1. Introduction 

In the large-scale Universe, one aspect of the search for regularities is connected by testing 
the radial velocity of galaxies. These velocities can be observed as having arbitrary values, 
regular patterns regarded as periodisation, discretisation or quantization of galaxy redshifts. 
The discretisation of redshift for astronomical objects can be discussed independently for three 
cases, namely galaxies, quasars and large-scale periodicity (120 Mpc). The latter studies have 
been discussed (Bajan et al. 2003), so we will not repeat these investigations here. In our 
previous paper (Bajan et al. 2004), the quest for quasar redshift periodicity is described. Here, 
only the main points from its history, together with some recent results, will be mentioned. 

The story of redshift periodisation started when Burbidge (1968) noted the existence of 
sharp peaks in redshift distribution, grouped around the values of z = 0.01 and z = 1.95. 
He also found periodicity in the redshift distribution, which can be described by the formula 
z = 0.01-n. 

The existence of periodicity in the log(l + z) scale was established by Karlsson (1971). 
In early stage of these investigations, association between bright low-redshift galaxies and 
quasars was sought in order to check for the non-cosmological origin of QSO redshifts, as well 
as to test for the possibility of ejection of QSO from parent galaxies, a problem investigated 
even now (Bell 2004). The link between GRBs and massive star formation is also suggested. 
Over the next years, the effect of periodisation was either confirmed or denied on the basis of 
incorrectly applied statistics and/or possible selection effect. The selection effects have been 
discussed by many authors. Some of them, (e.g. Karlsson 1971), claimed that the observed 
redshift distribution is not due to the selection effects, while other authors (e.g. Basu 2005) 
expressed the opposite opinion. 

Burbidge and Napier (2001) took into account a new sample of high redshift quasars lo- 
cated in projection on celestial spheare close (< lO") to the low redshift galaxies. The exis- 
tence of periodicity in the log(l -|- z) scale, was confirmed. The 2dF QSO Redshift Survey 
containing over 10000 objects and the 2dF Galaxy Redshift Survey with over 100000 galaxies 
served as an objective basis in the periodicity search for quasar-galaxy pairs. Having investi- 
gated 1647 objects, Hawkins, Maddox & Merrifield (2002) found no periodicity. Bell (2004) 
used two quasar samples, one with high redshifts of 2.4 — 4.8 and the other with low red- 
shifts of 0.02 — 0.2. He showed that all peaks in these two redshift distributions occur at the 
previously predicted preferred values. 

Napier & Burbidge (2001) re-examined the redshift distribution in the 2dF QSO Redshift 
Survey, detecting periodicity. Arp, Roscoe & Fulton (2005) showed that the maxima in redshift 
distribution in the 2dF and DSS surveys fit well into the formula. Bell's (2004) investigation 
showed clear periodicity in the redshift distribution of QSO. The latest work on redshift peri- 
odicity was done by Basu (2005). He investigated the same 33 objects (GRBs, QSO and active 
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galaxy) as Burbidge (2003) and contrary to Burbidge, he claimed that all existing peaks are 
due to observational and analytical selection effects. 

The discussion of galaxy redshift quantization started with the work of Tifft (1976). He 
claimed that the redshifts of galaxies in the Coma cluster were preferentially offset from each 
other in multiples of 72.46 km • s~^. A few years later, the existence of global periodicity was 
reported (Tifft & Cocke 1984); however, the period was not 72km • s~^, but 36km • s~^ or 
possibly 24km • s~^. It should also be noted that the values of 24km • s~^, as reported by 
Tifft & Cocke (1984), were not confirmed by the later investigations. 

We know that this subject is not popular and usually very suspicious at first glance. How- 
ever, on the basis of the claimed results of redshift periodisation, several theoretical papers 
pointed out the necessity for so-called new physics were published. 

Therefore, we decided to check if the discretisation does occur. We share the opinion 
expressed by Hawkins et al. (2002) that all these effects should be carefully checked. They 
claimed: 

"The criticism usually leveled at this kind of study is that the samples of redshifts have 
tended to be rather small and selected in a heterogeneous manner, which makes it hard to assess 
their significance. The more cynical critics also point out that the results tend to come from a 
relatively small group of astronomers who have a strong prejudice in favour of detecting such 
unconventional phenomena. This small group of astronomers, not unreasonably, responds by 
pointing out that adherents to the conventional cosmological paradigm have at least as strong 
a prejudice towards denying such results. 

We have attempted to carry out this analysis without prejudice. Indeed, we would have 
been happy with either outcome: if the periodicity were detected, then there would be some 
fascinating new astrophysics for us to explore; if it were not detected, then we would have 
the reassurance that our existing work on redshift surveys, etc, has not been based on false 
premises." 

Iwanowska (1989) claimed that the spatial distribution of galaxies belonging to the Local 
Group and globular clusters located close to our Galaxy is linear, i.e., these objects form long 
chains. For galaxies located in such lines, Zabierowski ( 1 993) distinguished 5 different redshift 
groups finding that the observed periodicity is consistent with the Tifft' s (1976) value of 72 km- 
s~^. He counted mean velocities in each group, using old data (RC3 (de Vaucouleurs et al. 
1991). Rudnicki, Godlowski & Magdziarz (2001) addressed the same problem, using better 
observational data. They considered 40 galaxies as well as globular clusters and performed a 
simple statistical analysis based on the calculation of the mean values of redshifts in bins and 
their dispersions. They tested strict quantization, that is precise multiplication of the value of 
36 km • s~^, finding no effect. However, with the result of the power spectrum analysis, the 
weak effect of periodisation, i.e. the grouping around some values of galaxy radial velocities, 
was noted for galaxies situated in two of the Iwanowska's lines. 
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The expected interpretation of galaxy redshift quantisation is that redshifts are not cos- 
mological. Moreover, there are suggestions that clusters of quasars evolved into clusters of 
galaxies, so the galaxy redshift distribution should reflect this fact. 

Tifft (1996) claimed that galaxy redshift distribution is quantised, because time is quan- 
tised. If this effect global, it means that it should be observed in aU galaxy structures. 

The number of galaxies in the Local Group is small but in the present work we consider 
aU galaxies regarded as members of the Local Group (hereafter referred to as the LG) (Irwin, 
2000; van den Bergh 1999). We apphed the power spectrum analysis as the statistical tool in 
our investigation. 

We used the standard method, i.e. power spectrum analysis, corrected for its possible dis- 
advantages when apphed to this problem. Moreover, we adopted a totaUy new approach. We 
use the Harm function as a weighting factor together with the jackknife error estimation, first 
used by Hawkins et al. (2002) (who investigated quasars, not galaxies). Because it is quite a 
new approach, we are enclosing several plots showing not only the results but also detailed 
analysing properties of the method. Such detailed analysis was not performed in the Hawkins 
et al. (2002) paper. Our result is that a weak effect of periodisation is observed for galactocen- 
tric velocities, while for heliocentric and LG-centric radial velocities no effect was observed. 

The paper is organized in the following manner. The second section presents observational 
data, and the third section describes the method of analysis applied. The fourth section consists 
of our results. In the last section we give our conclusions. 

2. Observational Data 

There is no common agreement as to which galaxies belong to the dynamical aggregate caUed 
the Local Group (LG). We considered 55 objects (see Tables 1 and 2) in our vicinity taken from 
the Irwin's list (Irwin 2000) (based on Mateo (1998)) together with 7 galaxies, mostly within 
the Maffei group, which could also probably be regarded as the LG members (Iwanowska, 
1989) (marked in Tables 1 and 2 as M). Considering various parameters, van den Bergh (1999) 
concluded that only 32 objects can be LG members, while 3 further objects can be regarded 
as possible LG members. We decided to perform aU calculations using these two sets of data. 
Sets based on the Irwin List were denoted as A, while those based on the van den Bergh's 
(1999) were denoted as B. 

It should be noted that in the van den Bergh hst there are 7 galaxies (not including Phoenix) 
without redshift, while the Irwin's hst contains 2 more such objects (Cetus and Cam A). In this 
manner, 46 and 28 galaxies remain to be analysed in each sample respectively. We denoted 
these samples as I and II. Separately, we analysed pure Irwin data (39 galaxies) as sample III. 

The search for any systematic effects in the distribution of redshifts requires precise knowl- 
edge of redshifts. There are some discrepancies among various redshift determinations. In or- 
der to avoid the influence of this factor on the results, we analysed separately samples with 
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radial velocities taken from the Irwin list (2000) (sample I) from those from van den Bergh's 
(1999) (sample II). Whenever no data was available for sample I, we decided to use the data 
from RC3 (de Vaucouleurs et al. 1991). 

The latest data allows us to find redshifts for 8 of the total 9 galaxies (marked as N in 
Tables 1 and 2), these redshifts were not known prior to these analysis. Redshifts for 7 of 
them are noted by Karachentsev et al. (2002). We repeated the analysis with these objects 
considered. However, we decided to exclude Tucana because its radial velocity measurement 
is uncertain. It should be noted that for the 6 remaining galaxies redshifts are also included in 
the new version of the Irwin's list (sample IV). Moreover, for two galaxies, the new version 
of the Irwin's Hst replaces the old redshift(NGC 147 ISTkms^^ andNGC 221 190 kms^^ 
). In our analysis we do not include. Cam A either (noted in NED as uncertain, with no errors 
determined). Finally, adding 6 galaxies to the Irwin's and the van den Bergh's lists we obtain 
samples of 45 galaxies (sample IV) and 34 galaxies (sample V) respectively. A detailed list 
of galaxies is given in Tables 1 and 2. Column (10) gives heliocentric radial velocity for set 
A (samples I, III and IV), while column (11) for set B (samples II and V). As seen from the 
data, errors of measuring radial velocities are really small for the majority of galaxies. 

The heUocentric radial velocities of galaxies should be corrected for the motion of the Sun 
relative to the center of our Galaxy and/or LG. This is usually done by applying the correc- 
tion to the center of our Galaxy only. There are several prescriptions for how to perform this 
reduction. In this paper, we applied various galactocentric corrections known from literature. 
We analysed the following solar motions: 

a): velocity v = 232k m • s^^ in the direction of / = 88°, b = 2° as proposed by 



Guthrie & Napien ( 1199 11) denoted as a, 

hy.OTv = 233±7krn-s~^ in the direction of / = 93° ± 10°, 5 = 2° ± 5° also proposed by 



Clube & Waddingtonld 19891) denoted as 5, 

c): and = 213 ± 10 k m • s^^ in the direction of ^ = 93° ± 3°, 6 = 2° ± 5° according to 



Guthrie & Napien ( 1 19961) denoted as c. 



There was also a suggestion that the local standard of rest has an additional radial com- 
ponent directed outward our Galaxy (Clube & Waddington 1989). We decided to add this 
component into correction a. In such a way, we obtained a new correction denoted as d with 
the following values: v = 234 km • s-\ I = 98°, 6 = 2° (Guthrie & Napier 1991). Pure 
heliocentric data that is without any corrections hereafter denoted as e were analysed, too. 

Moreover, the correction of the Sun's motion with respect to the LG centroid was included 
(Coutreau & van den Bergh 1999). This gave the velocity of the motion of the Sun relative 
to the LG of u = 306 ± 18km • s^^ toward an apex at ^ = 99° ± 5° and b = -3.5° ± 4° 
(/). Two additional corrections for the motion of the Sun relative to the LG were also taken 
into account. The correction denoted as g with the value of f = 305 ± 136km • s^^ toward 
I = M° ± 48°, b = -34° ± 29 calculated (Rauzy & Gurzadyan 1998) from two subgroups 
(ours and that of M31) in the LG, and "the historical" velocity (Yahil, Tammann & Sandage 
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1977), denoted ash:v = 308 ± 23 km • s'^ directed toward I = 105° ± 5 and 6 = -7° ± 4'' 



It is interesting to note that the spatial structure of the Galaxy is flat. Therefore, the correc- 
tion to galactocentric radial velocities hes in the Galactic plane. Therefore, we checked how 
the situation changes when a fictitious perpendicular component of the Sun velocity vector is 
assumed. We considered a velocity of w = 224 km • s ~ ^ directed toward I = 109° and 6 = 65° 
as correction i. FinaUy, the new proposed correction of the Sun's motion (Karachentsev et al. 
2002) in respect to the LG of w = 316km • s~^ toward the apex at Z = 93° and 6 = -4°, 
denoted as j, was used. AU of the corrections are summarized in Table 3. 

3. Method of analysis 

The power spectrum analysis (PSA) (Yu & Peebles 1969; Webster 1976; Guthrie & Napier 
1990), together with the Rayleigh test (Mardia 1972; Batschelet 1981), has been used as a 
statistical tool. It was shown (Newman, Haynes & Terzian 1994) that this method is very useful 
for finding periodicity among irregularly distributed points. The Rayleigh's test (Mardia, 1972; 
Batschelet 1981) is a simple test of uniformity which aUows one to detect periodicities in 
irregularly distributed points. For a given frequency, the Rayleigh power spectrum corresponds 
to the Fourier power spectrum, as well as measuring the probabiUty of the existence of a 
sinusoidal component. 

Let us assume that m points are distributed along a finite fine with coordinates Xi, where 
i= 1, . . . , m. We can define the phase with respect to the period P: 



The phase which corresponds to the i-th x coordinate unambiguously describes the radius 
vector 7j. 

The length of the vector R, constituting the sum of all li vectors, can be used for testing 
isotropy of the distribution of points Xi with respect to the period P 



From the formal point of view, the vector R describes a random walk in the plane after m 
steps. 

We use the statistic R^ defined by: 



were used. 



$i = 27rXi/P. 



(1) 



JV 




(2) 




(3) 



The probabiUty distribution of R"^ 



can be calculated from the null hypothesis: 



$G (0,27r), 



(4) 
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where is the probabiUty that phase corresponding to point Xi is located in the 

interval ($, $ + The distribution of corresponds to the Fourier power spectrum for 
function /($) = - 
It is known that the variable: 



where Jo is a Bessel function. This formula is also valid for a small number of objects, which 
is the case of the LG. 

The distribution p{s, m) is calculated numerically by integrating approximations of the 
Bessel function using the Romberg method (Press et al. 1992). For large m, it could be also 
aproximated as a distribution with 2 degrees of freedom. 

Error bars of s{P) can be estimated using the "jackknife" technique of drawing all possible 
samples oiN — 1 values from N data points, repeating the power spectrum analysis on these 
samplings. Such a procedure allowed us to calculate the standard deviation in the derived 
values of s crj(P). The best estimator for the standard errors in the value of s is then just 
s/N - 1(7 j (Hawkins et al. 2002). 

The simulation of the power spectrum for random uniformly distributed data is presented 
in Figure 1 . The diagrams showing the values of s-statistics versus 1/P present the results of 
the power spectrum analysis. There are several peaks in each diagram. These peaks allow one 
to find each possible period, as well as to investigate the significance of each particular peak 
in the power spectrum. The level of significance of each peak is given by C = 1 — (1 — p)"', 
where p is the probability of obtaining, from the theoretical, random distribution the value 
of the s-statistics equal to or greater than the observed value of the s-statistics, while rit is 
the number of independent peaks within the analysed frequency range (Lake & Roeder 1972; 
Guthrie & Napier 1991). 

The additional test increasing the efficiency of the test for weak clustering, following 
(Webster 1976; Scott 1991), is based on the summation over the whole power spectrum. This 
sum gives the value of SI = ^Si, with a distribution with 2nt degrees of freedom. Thus, 
the expected value of the SI statistic is 2nt. The S'Z-test can be used for testing ihe random- 
ness of the distribution. 

The clustering statistics Q is equal to the value of the SZ-statistics over its expected value. 
The expected value for Q statistics is calculated in the case of a random walk (2nt). For a 
random distribution, the expected value of Q is equal to 1 , with the error of : a (Q) = 1/ y/n^ 
(Webster 1976). We tested the hypothesis that the value of Q is greater than unity rather than 
equal to this value. In our analysis, we decided to consider the first 50 peaks. In this case, at 
the significance level of 0.05, the critical value of the SI statistics is 124.3, with 135.8 for the 
significance level of 0.01, also note that a (Q) = 0.14. 



s(P) = 2R^/m, 



(5) 



has the following distribution (Webster 1976): 
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Newman et al. (1994) pointed out that Yu and Peebles' version of PSA can be applied 
correctly only when a uniform distribution function is tested. As seen from Figure 2, such an 
approximation could be accepted for raw data but not when a correction for the solar motion 
is included. 

Hawkins et al. (2002) discussed the power spectrum method when the data is not uniformly 
distributed. They proposed to use the window function and showed that the power spectrum 
method works well in that case. Now, the power of s at the period of P is given via the 
formulae (Hawkins et al. 2002): 

s{P)=2B?/Y,wl (7) 

where 

N N 

R^{P) = C^wicos^i)^ + C^WiSin^if, (8) 

i=l i=l 

Following Hawkins et al. (2002), we repeat our analysis using the Hann's function as a 
weight: 

/ 27rxA 

where L is chosen to cover the whole range over which the data is selected. 

It should be pointed that now the expected value of the clustering statistics Q is not nec- 
essarily equal to 1, especially for a small number of points. Therefore, we decided to test 
how the PSA with the Hann's weighting function is working. The observed distribution of 
galactocentric redshift in the LG is normal (Fig. 2). Using the K-S tests, we checked that the 
distributions denoted as 1 and 11 and containing 46 and 28 galaxies respectively, are normal at 
the significance level of a = 0.05. Thus, in the future analysis we use the Gaussian redshift 
distribution, as theoretically expected distribution. 

We run 1000 simulations for data distributed normally, and with variance taken from the 
real data with correction a. We find that in the case of sample 1 (46 galaxies), the mean value of 
Q is equal to 1, but = 0.195. For sample 11 (28 galaxies), the mean value of Q decreases 
and is equal to 0.90, but a{Q) = 0.177. This means that the noise value of s = 2 (see Hawkins 
et al. 2002) is slightly changed. The example of the simulation of the power spectrum for 
normally distributed data weighted using the Hann's function is presented in Figure 3. One 
can see that when the data are apodized, the expected height of the peak decreases. 

4. Result 

We started our investigations with a classical power spectrum analysis (Yu & Peebles 1969; 
Webster 1976). The two basic samples of galaxies (based on sets A and B), containing 46 
and 28 objects each, denoted as sample I and n respectively, were analysed. For each sample, 
10 different reductions of heliocentric radial velocities to galactocentric velocities or a LG 
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centroid, were performed. In Figure 4 the power spectra of sample I analysed as s versus 
1/P diagrams in the range of 1/P G (0.0025, 0.05) is shown. The visual impression is that 
some peaks are clearly distinguishable, which can be regarded as concentration around some 
particular values of the period. These values are close to 36 km- s~^, and 24 km- s~^. However, 
this impression is not confirmed by statistical analysis. 

The result of our analysis is given in Table 4 (sample I) (and Table 5 for sample II respec- 
tively). The first column of both tables describes the analysed sample, while the next two give 
the values of the velocity v at which the most prominent peak was observed, together with 
the probability P{s) that the peak is generated by a random distribution (when we restrict for 
period P G (20, 100) km - s'^). 

From Table 4, it is easy to see that the most prominent peaks are statistically not significant, 
at the significance level of 95%. Therefore, we concluded that in random distribution such 
peaks can also appear. The power spectrum analysis shows that the peaks observed in the 
period distributions are consistent with the assumption of randomness. 

The second aim of the power spectrum analysis was to check for the possibility of weak 
clustering using SI and Q-statistics. The result of our analysis is also given in Table 4. The last 
two columns of the table contain the values of SI and Q statistics. It can be easily seen that in 
the case of galactocentric correction prescriptions a, h, c, both 5*/ and Q-statistics confirm the 
hypothesis that the redshift distributions are non-random at the confidence level of 95%. For 
prescriptions a and b this seems to be confirmed even at the level of 99%. 

For the prescription d computed using the probable velocity of the local rest frame, we do 
not find any departure from isotropy. A comparison with the previous cases allows us to state 
that this vector does not correspond to reality. The reduction to the LG centroid, using any of 
our vectors, did not reveal any statistically significant individual preferred values of velocities 
and showed that the distributions are random. 

The high value of s given by the PSA is expected at scales greater than the clustering scale 
(Yu & Peebles 1969; Scott 1991). We deal with one structure (Local Group) only, so this is 
not the cause of the reported non-randomness. Nevertheless, we noted that from the theoretical 
point of view the non-random distribution detection by the PSA can be due to the inclusion of 
some external objects. In order to examine this possibiUty, we performed three tests: 

i) we consider sample II, from which all objects with untypical values of redshift were 
excluded (van den Bergh 1999), 

ii) we repeat our analysis of sample I without the Argo galaxy, which is its most pro- 
nounced external member, 

iii) because the artificially high values of s given by the PSA are observed mainly in 
the first mode (Webster 1976), for the period we perform the PSA also in the interval i = 
2, . . . , 51 instead of the first 50 peaks (this is also the reason why we restrict our analysis to 

Pg (20, 100)km-s-i). 
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All these changes alter the statistics values, but for the prescriptions a and b, the PSA 
distribution remains non-random. 

Karachentsev et al. (2002) showed that, for galaxies belonging to the Local Group, errors 
in heliocentric velocities are generally ±5 km • s~^ or smaller. This is true for galaxies belong- 
ing to sample n. However, it should be noted that some of the galaxies belonging to sample I 
have large errors in their reported radial velocities. This is the case for the sample of 7 galax- 
ies (marked in Tables 1 and 2 as M) which Iwanowska (1989) regarded as the probable LG 
members (although this is still very questionable). It should be pointed out that large errors in 
the data can mask any existing periodisation (however, this cannot be a source of a fictitious 
signal in PSA). Now we repeated our analysis for the pure Irwin's list, withouth those 7 galax- 
ies (sample III). The result (Table 5) is similar to that obtained for sample I but the value of 
statistics SI for reduction a, b and c is even larger than in sample I. This confirms that these 7 
galaxies denoted as M in Tables 1 and 2 should be rejected from further analysis. 

When we add 6 new galaxies (marked in Tables 1 and 2 as N), we obtain (for samples IV 
and V analysed now) that the PSA distribution clearly remains non-random for prescriptions 
a only; however a weak effect survives for prescriptions b and c as well (see Table 5). 

Still, to be sure that the obtained effect is real, we should take into account the fact that 
our data are not uniformly distributed. In that case (Hawkins et al. 2002), the PSA must take 
into account the data window function. Therefore, we repeated our analysis using the Hann's 
function as a weighting function (Hawkins et al. 2002). For sample II, we clearly obtain that 
the distribution of the PSA remains non-random for reduction a (Table 6, statistic Q). Figure 
5 also shows that, for prescriptions c, the most prominent peak is significant at about 2a level 
above the noise value. For prescriptions a and b, the significance of the most prominent peak 
is reduced below the 2(j level. A similar situation is seen in Figure 5 for sample V, where for 
prescriptions a and c the peak is significant at about 2a level. It should be noted that the PSA 
with an appropriate window function does not show any signs of periodicity in other cases. 
This suggests significant periodisation close to 24 km • s~ ^ . 

Moreover, it should be pointed out that using the window function we change not only 
the expected value of the clustering statistics Q, but also the expected height of the peak. 
As a result, the probabihty P{s) formally obtained on the basis of Equation 6 (presented in 
Table 6) is not vaUd any longer. Again we run 1000 simulations for the normally distributed 
data with variance taken from the real data. The histogram for the distribution of the height 
of the most significant peak with data weighted using the Hann's function for sample 11 (28 
galaxies) with variance taken from the real data with correction a, is presented in the left panel 
of Figure 6, the cumulated distribution is presented in the right panel. We can see that peaks 
with 2B? /m > 8.2 are significant at the 95% level. Please note that, at least for prescription 
c, for samples 11 and V the most pronounced peaks are significant at the level of 95%. A 
similar situation occurs for sample V, prescription a. This confirms weak periodisation close 
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to 24 km • s and is in contrast with the previous cases (without using the data window 
function), where all peaks were not significant. 

5. Conclusions 

In our analysis, we considered the nearby galaxies which are widely accepted to be LG mem- 
bers. Moreover, we defined a sample of galaxies according to several physical properties and 
not only to spatial distribution, as belonging to the LG. We took into account several possible 
ways of correcting the hehocentric velocity to the galactocentric velocity. The correction to 
the LG centroid was also considered. We used the best data currently available. 

The preliminary statistical analysis excluded the possibility of strict redshift quantization 
in the LG (Rudnicki et al. 2001). The lack of strict multiplicity of the value of 36 km • s~^ 
based on good data does not confirm Zabierowski's claim (Zabierowski 1993) to this effect. 

We obtained that distributions of galaxy redshift seem to be non-random when the correc- 
tion for the motion of the Sun relative to the center of our Galaxy is taken into account. The 
argument revealing a possible significance of the result is the fact that it is practically inde- 
pendent of the samples. There is some periodisation though which in some cases is close to 
36 km • s~^, while in other cases it is close to 24 km • s~^. However, this effect is statistically 
weak. The statistical significance of our result as seen from the Tables is rather low. 

The most important argument is that the result survives the PSA when the data window 
function is taken into account. This provides that the effect found is real. However, only pe- 
riodisation close to 24km • s~^ seems to be confirmed at the significance level of 95%. One 
can see that the strongest effect is obtained for galaxies from the van den Bergh's list. Van 
den Bergh (1999) concluded that only 32 objects can be LG members, while 3 further objects 
can be regarded as possible LG members. The Irwin's Ust contains more galaxies for whom 
membership in the LG is questionable. 

With correction for the Sun's motion with respect to the LG centroid, we obtained no 
periodisation. We did not obtain any statistically significant preferred values of velocities in 
this case. It is interesting to note that the spatial structure of the Galaxy is flat. Therefore, the 
correction to the galactocentric radial velocities lies in the galactic plane. We also checked 
how the situation changes when a fictitious perpendicular component of the Sun's velocity 
vector is assumed. In this case there is no periodisation. This result, together with the analysis 
with prescription e (pure heliocentric velocity) shows that the effect is diluted for an artificial 
value of this vector. Thus, this provides evidence that the observed weak periodisation is not 
an artifact of the galactocentric correction. 

The latest data allows us to include redshifts for 6 galaxies determined by Karachentsev et 
al (2002). When we repeated the analysis with these objects, the effect decreased. However, it 
should be pointed out that these objects are cormected with M31 (and they are probably M31 
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satelites). Thus, any possible periodisation could be stronger than for other galaxies masked 
by the local influence. 

The interpretation of this result is neither unique nor clear. We would like to point at the 
following possibilities: 

1 . The effect of periodisation of galaxies in the LG can be connected artificially with the 
manner of galactocentric correction due to either velocity determination or to the data itself 
but the above-mentioned diminishing of non-randomness when artificial velocities are used 
does not seem to support this conclusion. 

2. Periodisation is found only when galactocentric velocities are included, and it is not 
connected with the LG centroid. This is possible, for example, when the effect is global as 
claimed by Tifft and it is not connected with the LG. 

3. Periodisation is a relic of quantum effects in the early Universe (Tifft 1996). Thus, the 
effect is not connected with the LG, which was formed during much later epoch. 

4. Periodisation is a real phenomenon, also connected with the LG, but the correct LG 
structure is not known at present. 

Clearly, power spectrum analysis is a good tool for studying non-randomness in velocity 
data. The main result of the paper is that the distributions of galaxy redshift seem to be non- 
random, but each specific individual particular value of periodisation is statistically marginal. 
Only periodisation close to 24 km • s^^ seems to be confirmed, and only at the significance 
level of 95%. One of the reasons is that errors in the data are large. However, it should be noted 
that in the version of the PSA considered, non-systematic errors can not produce "false peri- 
odisation" but only "destroy" any real existing effect. Measurement errors will be important if 
explicitly taken into account. In such a version of PSA, we analysed the sample from the point 
of view of the objects likelihood rather than using its best values. However, such a procedure 
changes the initial assumptions of the method as well as the analysed statistics, therefore it is 
not considered in our present paper It will be analysed in our future work. Further investiga- 
tion should concern such large-scale structures as the Local Supercluster, Coma/A 1367, the 
Perseus or Hercules Superclusters, when more accurate HI velocity data become available. 
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Fig. 1. Example of simulation of the raw 
power spectrum of a random uniform distri- 
bution — sample II (28 galaxies). Dashed 
lines show errors derived by applying the 
jackknife estimator. 



Table 3. Motion of the Sun. 
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Fig. 3. Example of simulation of the power 
spectrum for normal distribution derived with 
data weighted using the Hann's function — 
sample II (28 galaxies). Dashed hnes show 
errors derived by applying the jackknife es- 
timator. 



Table 4. Statistical analysis of samples inves- 
tigated in the Local Group (Sample I) 
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Fig. 2. Redshift distribution histogram for raw data (left side), and those corrected for the solar 
motion (correction a). 
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Fig. 4. Raw power spectrum of sample 1 (46 galaxies). Dashed lines show errors derived by 



applying the jackknife estimator. 
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Fig. 5. Power spectrum of sample 11 (28 galaxies), III (39 galaxies) and V (34 galaxies) after 
apodization with the Harm's function . 
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Fig. 6. Histogram for distribution of the height of the most significant peak with data weighted 
using the Hann's function — sample II (28 galaxies) (left panel), and cumulated distribution 
of the height of the most significant peak with data weighted using the Hann's function — 
sample II (28 galaxies) (right panel). 
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Table 1. List of galaxies. 
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Table 2. List of galaxies. 
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Table 5. Statistical analysis of samples investigated in the Local Group (Samples II - V) 
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Table 6. Statistical analysis of samples in- 
vestigated in the Local Group (with data 
weighted using the Hann's function). 
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